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Magnetoresistance oscillations in two-dimensional electron systems driven simultaneously by a 
strong direct current and a microwave irradiation, are analyzed within a unified microscopic scheme 
treating both excitations on an equal footing. The microwave-induced resistance oscillations are 
described by a parameter e w proportional to the radiation frequency, while the dc-induced resistance 
oscillations are governed by a parameter ej proportional to the current density. In the presence of 
both a microwave radiation and a strong dc, the combined parameter e^ + ej is shown to control the 
main resistance oscillations, in agreement with the recent measurement [Zhang et al. Phys. Rev. 
Lett. 98, 106804 (2007)]. 
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Interest in the direct current effect on magnetore- 
sistance in two-dimensional (2D) electron system (ES) 
has recently been revived due to the exciting experi- 
mental observations that in the case without irradia- 
tion, a relatively weak current can induce substantial 
magnetoresistance oscillations and negative differential 
resistivity*^!^ Recent careful measurements^ and theo- 
retical analysis^ disclosed accurate oscillating charateris- 
tic having the period Aej ~ 1, expressed by a dimension- 
less parameter ej proportional to the ratio of the current 
density to the strength of the magnetic field. 

The main oscillations appearing in the linear mi- 
crowave photoresistance were known to be periodic in 
e u = lu/uj c (lu is the angular frequency of the radiation 
and oj c = eB /m is the cyclotron frequency of an effective- 
mass m electron in the magnetic field B), having the pe- 

riod Ae w = 1^8^0,11,12,13,14,15,16,17,18,19,20 Recent stud _ 

ies by Zhang et alZ^ on magneto-photoresistance in a 
2DES subject to a microwave irradiation and a finite 
direct-current excitation revealed that the maxima (min- 
ima) of the differential resistance can evolve into minima 
(maxima) and back, and the oscillations are governed by 
a parameter e — e u + ej , combining the dc related param- 
eter ej and the radiation- frequency related parameter e u . 

The preliminary explanation of the observed resistance 
oscillations under simultaneous dc and ac excitations, 
which is based on indirect electron transitions consisting 
of a radiation-induced jump in energy and a dc-induced 
jump in space between Hall-field-tilted Landau levels^! 
is obviously not satisfied. A systematic theory capable of 
treating both excitations simultaneously within a single 
framework is highly desirable. 

It is better to treat the magnetotransport subject to a 
strong dc excitation with a scheme direct using the elec- 
tric current, rather than field, as the basic control pa- 
rameter. Such a theory has been developed) 15 > 19 In this 
theory the integrative drift motion of the electron sys- 
tem opens new channels for individual electron to transit 
between different Landau levels through impurity and 
phonon scatterings, thus the direct current and the mi- 
crowave radiation are dealt with on an equal footing. 



We consider a 2DES having N s electrons in a unit area 
of the x-y plane subject to a uniform magnetic field B — 
(0, 0, B) in the z direction. When an electromagnetic 
wave with incident electric field E- ls sinwi, irradiates the 
2D plane together with a dc electric field Eq inside, the 
steady transport state of the system can be described by 
the electron drift velocity Vq and an electron temperature 
T c , satisfying the force and energy balance equations; 19 

N s eE + N s e(v x B) + F = 0, (1) 
N s eE Q ■ v + S p - W = 0. (2) 

Here, the frictional force resisting electron drift motion, 

oo 

F = ^|L%||)| 2 q\\Jn(On 2 (q {l) u> + nw), (3) 

qy n— — oo 

is given by the electron density correlation function 
Il2(q\\,f2), the effective impurity potential U(qn), a 
radiation-related coupling parameter £ in the Bessel func- 
tion Jn(£)> and ujq = qn -Vg. The electron energy absorp- 
tion from the radiation field, S p , and the electron energy 
dissipation to the lattice, W, are given in Ref.[l^. The 
nonlinear longitudinal resistivity and differential resistiv- 
ity in the presence of a radiation field are direct obtained 
from Eq. |T]) by taking Vo (i.e the current J — N s evo) in 
the x direction, vq = (vo, 0, 0) and J — (J, 0, 0): 

R xx = -F /(N*e 2 v Q ), (4) 
r xx = = -{dF Q /dv Q )/{N 2 s e 2 ). (5) 

In the absence of microwave radiation, the formulation 
reduces to those presented in Ref.[fJ. The initial suppres- 
sion and subsequent oscillation of the magnetoresistivity 
arise from the current-assisted electron transitions and 
are controlled by the parameter 

_ uij 2mkpVo / 87T m J 
6j "~ c = eB = \IN~ S ~^B' 

The oscillating valley-peak pairs in the R xx -vs-ej curve 
appear around the positions ej = r\m [m = 1,2,...), with 
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r] ~ 1, depending on the form of the scattering potential 
and on the electron temperature. 

In the case of linear photoresistance, vq — > 0, the elec- 
tron transition (intra- and inter-Landau levels) can take 
place by absorbing or emitting n photons of frequency 
lu and jumping from Landau level I to I' (across to lev- 
els). This process gives rise to terms of ±n in the Mo- 
tional force Fq expression Eq. |(3J| . The periodicity of the 
electron density correlation function Il2(q»,0 + lo c ) = 
//2(<7|| , (2) at low electron temperature and many Landau- 
level occupation, leads to the appearance of peak-valley 
pairs in the linear magnetoresistivity around the posi- 
tions of nuj = moj c £Q The primary peak- valley pairs 
come from the single-photon process, n = 1, showing 
up around the node positions to — moj c , or 

e u =co/w c = m= 1,2,3, (7) 

and the primary period of R xx oscillation is Ae w = 1. 
Two- and multiple-photon processes (n > 2) yield sec- 
ondary peak-valley pairs^ 

When there is a finite bias current in 2DES, in view 
of the extra energy c^o = Q\\ • vq provided by the elec- 
tron drift motion together with the energy nw provided 
by the absorption or emission of n photons of the radi- 
ation field, an electron can be scattered by elastic im- 
purities and jumps across to Landau levels of different 
energies. The resonance condition in this case should 
be determined by r)U)j + ncu = rriLu c . Focusing only on 
single-photon processes (n — 1), the valley-peak pairs of 
nonlinear magnetoresistance oscillations are anticipated 
to appear around the positions 

^w + V e j = e u> (1 + VUj/u) = m (8) 

or 

e w = uj/uj c = m/(l + 7j), (9) 
where m — 1, 2, 3, • • and 

Jj=T}Dj/u. (10) 

Therefore, the oscillating behavior of R xx is controlled 
by the parameter e u + r)6j, or e w (l + 7,). 

Figure 1 presents the calculated longitudinal mag- 
netoresistivity R xx vs e u = lu/lu c for a GaAs-based 
2DES with carrier density N s = 3 x 10 15 /m 2 , 
low-temperature linear mobility /iq = 2000m 2 /Vs 
at lattice temperature T = IK, irradiated by a 
linearly rr-direction polarized microwave of fre- 
quency w/27r = 100 GHz and incident amplitude 
E- ls = 5V/cm, and subject to bias velocities 2vq/vf = 
0, 0.001, 0.002, 0.003, 0.004, 0.005, 0.007, 0.01, 0.05, and 
0.1 (vf is the Fermi velocity), corresponding to current 
densities J = 0.06,0.11,0.17,0.23,0.29,0.40,0.57,2.9, 
and 5.7 A/m or the bias parameter ojj/oj — 
0,0.05,0.1,0.16,0.21,0.26,0.36,0.52,2.6, and 5.2, 
respectively. The elastic scatterings are assumed due to 
short-range impurities and the Landau level broadening 
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FIG. 1: (Color online) Longitudinal magnetoresistivity R xx 
vs = u/uj c for a GaAs-based 2DES with iV s = 3 x 
10 15 /m 2 , ^0 = 2000m 2 /Vs and a = 10 at lattice tem- 
perature T = 1 K, under linearly polarized irradiation of 
frequency u/2n — 100 GHz and incident amplitude E[ s = 
5V/cm and subject to bias current densities of uJj/uo = 
0,0.05,0.1,0.16,0.21,0.26,0.36,0.52,2.6, and 5.2. The elas- 
tic scatterings are due to short-range impurities. The inset 
shows Rxx as a function of uJj /u) at three fixed magnetic field 
strengths, lu/uj c = 1.16, 1.80, and 2.22. 



parameter is taken to be a = 10J^ Linear (zero dc bias, 
LUj = 0) magnetoresistivity exhibits typical feature of 
radiation-induced magnetoresistance oscillations, with 
two-photon process slightly showing up as a shoulder 
around e w ~ 1.5. In the case of small dc bias, <C 1, the 
node positions of peak- valley pairs are shifted down from 
e u = m at zero dc bias to e u = m/(l + jj) w m — mjj. 
The shifted distances relative to the original positions, 
rn7j, are larger for larger to, but still locate within the 
same order of e w as long as mr/j < 1. When the shifted 
distance becomes as large as 0.5, the original peak will 
move to the range where previously there is a valley. 
In the case of strong dc bias, jj > 1, the original pairs 
of higher order (larger m) may shift to lower order 
range of e w . For instance, in the case of short-range 
impurity scattering 7] = 0.94, 6 for the bias current of 
loj/lo = 5.2 (7j = 4.9), the original pairs at e w = m 
with to = 1,2,3,4,5, and 6 are shifted to the positions 
e u = m/(l +7j) = 0.17,0.34,0.51,0.68,0.85, and 1.02. 
Except for the first one (e u — 0.17), which is beyond the 
range plotted, the other five pairs are clearly identified 
in the figure. 

The inset of Fig. 1 shows the change of the magne- 
toresistivity R xx with increasing bias current density J 
in terms of ujj/uj at three fixed magnetic field strengths 
u>/u> c = 1.16,1.80, and 2.22, respectively, near the first 
minimum, the second maximum and minimum of R xx . 
Here we show the fact that within a certain magnetic 
field range of resistivity minimum, R xx can be negative 
at small J but increases with increasing J and passes 
through zero at a finite Jr^ This is exactly what is re- 
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FIG. 2: (Color online) Longitudinal magnetoresistivity R xx 
(a) and differential magnetoresistivity r xx (b) vs e u = o;/o; c 
for a GaAs-based 2DES with iV s = 3.7 x 10 15 /m 2 and no = 
1200m 2 /Vs at T = 1.5 K, under irradiation of a linearly po- 
larized microwave of frequency lj/2tv = 69 GHz and incident 
amplitude E- 1S = 3.6 V/cm, and subject to bias current densi- 
ties of u)j/u) — 0, 0.19, 0.37, 0.56, and 0.93. The elastic scatter- 
ings are from a mixture of short-range (SR) and background 
(BG) impurities. 



quired for the instability of the time-independent small- 
current solution and for the development of a spatially 
nonunifor m 22 ! 23 or a time-dependent stated which ex- 
hibits measured zero resistance. Therefore, the region 
where an absolute negative dissipative magnetoresistance 
develops is identified as that of zero-resistance state. 

The calculated longitudinal magnetoresistivity R xx 
and differential magnetoresistivity r xx for another GaAs- 
based 2DES with carrier density N s = 3.7 x 10 15 /m 2 
and low-temperature linear mobility /io = 1200m 2 /Vs 
at lattice temperature T = 1.5 K are plotted in Fig. 2 
as functions of = u>/uj c , under the irradiation 
of a linearly x-direction polarized microwave of fre- 
quency oj/27t = 69 GHz with incident amplitude E- ls — 
3.6 V/cm, and subjected to bias drift velocities 2vq/vf = 
0, .002, .004, .006, and 0.01, corresponding to current den- 
sities J = 0,0.14,0.28,0.42, and 0.70 A/m or bias pa- 
rameter ujj/uj = 0,0.19,0.37,0.56, and 0.93, respectively. 
The elastic scatterings are assumed due to a mixture of 
short-range and background impurities and the broaden- 
ing parameter is taken to be a = 7. We see that po- 
sitions of peak-valley pairs and their maxima and min- 



ima of different e u orders m in the zero-bias R xx curves 
[Fig. 2(a)] are shifted downward with increasing current 
density u)j / uj following the same rule as stated in the case 
of Fig. 1. 

The oscillations of differential resistivity r xx exhibit 
shorter period and much larger amplitude than those of 
R xx , especially at high dc bias where R xx oscillations 
become less prominent, as shown in Fig. 2(b). Positions 
of peak-valley pairs and their maxima and minima in 
r xx are shifted downward further relative to those of the 
corresponding R xx at the same order m. Of particu- 
lar interest is that with increasing bias current density, 
in the e u range where the zero-bias magnetoresistivity 
(R xx = r xx ) exhibits positive maximum, the differential 
magnetoresistivity r xx can be driven down to a consider- 
able negative value, as indicated by the up-directed ar- 
rows in Fig. 2(b); and in the e w range where the zero-bias 
magnetoresistivity (R xx = r xx ) exhibits a minimum, the 
differential magnetoresistivity r xx can be driven up to 
exhibit a maximum, as indicated by the down-directed 
arrows in Fig. 2(b). 

According to electrodynamic analyses, a homogeneous 
state of current J is unstable if either the absolute dis- 
sipative resistivity R xx or the differential resistivity r xx 
becomes negativei 22 i 23 Thus, despite the detailed struc- 
ture yet to be explored, it is expected that, under appro- 
priate conditions, a zero-resistance state (ZRS) or some 
kind of quasi-ZRS could be measured in the region where 
R xx or r xx exhibits a negative value in a homogeneous 
microscopic analysis. With this in mind, the R xx and 
r xx behaviors in Fig. 2 could predict the following mea- 
sured results. In the vicinity around e u = 2.25, the ZRS 
showing up in the zero dc bias will disappear at bias 
u)j/u> > 0.10 because both R xx and r xx return to positive. 
Instead, one should observe a bulged r xx in the cases of 
ujj/uj = 0.19, 0.37, and 0.93 and a dented r xx in the case 
of tOj/uj — 0.56, as shown in Fig. 2. In the vicinity around 
e u = 1.25, where the negative R xx (thus the ZRS) main- 
tains up to the dc bias u)j/u> — 0.35, one may not detect 
r xx significantly different from zero at low dc biases but 
the bulged r xx should be observed at higher dc biases 
when R xx becomes positive and the ZRS disappears, as 
seen in Fig. 2 in the cases of ujj/u = 0.37 and 0.56. In the 
vicinity around e u — 1.75 and 0.8, the zero-bias magne- 
toresistivity R xx = r xx exhibits main peaks while the dif- 
ferential resistivity r xx becomes negative, respectively, at 
finite dc biases ujj/ui = 0.37 and 0.56, suggesting the pos- 
sible ZRS induced by the direct current. When further 
increasing current density, it disappears and becomes a 
bulged differential resistivity. In the range 0.6 < e u < 1 
the differential resistivity r xx at all three strong dc bi- 
ases (ujj/uj — 0.37,0.56, and 0.93) can become negative, 
implying possible resistance suppression or quasi-ZRSs 
induced by the bias current. These results indeed show 
up in a recent experimental observation.— 
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